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Directional Radiation and Photodissociation Regions 
in Molecular Hydrogen Clouds 

S. A. Balashev*, D. A. Varshalovich** , and A. V. Ivanchik*** 

Ioffe Physical- Technical Institute, Politekhnicheskaya 26, St. Petersburg, 194021 Russia 

Abstract 

Some astrophysical observations of molecular hydrogen point to a broadening of the velocity distribution for molecules 
at excited rotational levels. This effect is observed in both Galactic and high redshift clouds. Analysis of H2, HD, and 
CI absorption lines has revealed the broadening effect in the absorption system of QSO 1232+082 (z a b s = 2.33771). 
We analyze line broadening mechanisms by considering in detail the transfer of ultraviolet radiation (in the resonance 
lines of the Lyman and Werner H2 molecular bands) for various velocity distributions at excited rotational levels. The 
mechanism we suggest includes the saturation of the lines that populate excited rotational levels (radiative pumping) 
and manifests itself most clearly in the case of directional radiation in the medium. Based on the calculated structure 
of a molecular hydrogen cloud in rotational level populations, we have considered an additional mechanism that takes 
into account the presence of a photodissociation region. Note that disregarding the broadening effects we investigated 
can lead to a significant systematic error when the data are processed. 

Key words, high-redshift molecular clouds - absorption lines in quasar spectra - radiative transfer - QSO 1232+082. 



INTRODUCTION 



Molecular hydrogen, the most abundant molecule in 
the interstellar medium, is observed not only in our 
Galaxy QSavage et al. 1977| IShull fc Beckwith 19821 
Rachford et al. 2002) and Local Group galax- 
ies (jTumlinson et al. 20021 IBluhm et al. 2003|) but 
also in the spectra of high-redshift quasars (see 
Levshakov & Varshalovich 1985; Noterdacme et al. 2008). 



The importance of molecular hydrogen, particularly at 
the early evolutionary phases of the Universe, lies in the 
fact that it plays a key role in the theory of star formation, 
being the main coolant. Observations of various rotational 
H2 states allow the main physical parameters (kinetic tem- 
perature, pressure, particle density, etc.) to be determined 
in the molecular hydrogen clouds that existed 12 Gyr ago. 
The relative population of rotational H2 levels also makes 
it possible to estimate the intensity of ultraviolet radia- 
tion, whose high value can serve as an indicator of active 
star formation in early protogalaxies. In addition, preci- 
sion measurements of molecular hydrogen wavelengths at 
high redshifts are used to estimate the possible cosmolog- 
ical change in the proton-to-electron mass ratio (see, e.g., 
Ilvanchik et al. 2005|) . 

The ground X*E+ and excited B 1 ^ and C 1 ^ elec- 
tronic states are distinguished in the molecular hydrogen 
level structure. There are vibrational and finer rotational 
level structures in each electronic state (see Fig. [IJ. 

Since the dipole transitions are forbidden in the ground 
electronic state, molecular hydrogen is commonly observed 
using the system of lines between the ground X state and 
the excited B and C states, which are called the Lyman and 
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Werner bands, respectively. These bands lie in the ultravi- 
olet spectral range. However, for hydrogen at high redshifts 
(z ~ 2+4), the absorption lines fall within the visible spec- 
tral range. Therefore, it can be observed using groundbased 
telescopes with a high spectral resolution. 

Another peculiarity of the hydrogen molecule is the 
mechanism of its dissociation. The point is that the dissoci- 
ation continuum of molecular hydrogen is higher than that 
of atomic one. This means that there is no direct photodis- 
sociation of molecular hydrogen in the interstellar medium, 
because atomic hydrogen depletes the spectrum above 13.6 
eV. The dissociation proceeds through the mechanism that 
was first described bv lStecher fc Wil liams 19671 After their 
excitation in the Lyman or Werner bands, about 87% of 
the molecules return rapidly, in a time r ~ 10 -6 s, to the 
ground electronic state at upper ro- vibrational levels, while 
13% of the molecules fall into the continuum, i.e., dis- 
sociate ( Abgrall et al. 1992| . The undissociated molecules 
fall into the lower vibrational state through cascade tran- 
sitions and populate the long-lived rotational levels (see 
Dalgarno & Black 1974). This mechanism is called radia- 
tive pumping. 

The absence of dipole transitions in the lower 
vibrational, ground electronic state not only makes 
the rotational levels basically metastable but also 
strongly suppresses the direct radiative population. 
The typical physical conditions in molecular hydro- 
gen clouds are T ~ 50 + 150 K, n ~ 10 + 1000 sm~ 3 
(|Rachford et al. 20021 IShull et al. 2000|) . At such parame- 
ters, the collision rate is too low (compared to the rates of 
spontaneous transitions) for the excited rotational levels to 
be populated. Therefore, being in background ultraviolet 
radiation with an intensity I ~ 10~ 18 erg/sm 2 /s/Hz (for 
our Galaxy, at wavelengths of the order of A = 1000 A, 
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Figure 1. Level diagram for the hydrogen molecule. The ground X : E+ and excited B X E+ and C 1 !^ electronic states are 
presented. The vibrational structure of the electronic levels is shown. The rotational structure for the lower vibrational, 
ground electronic state, in which the times of spontaneous transitions are indicated, is presented on the right panel. 



IDraine 1978[) . the excited rotational J > 2 levels are 
populated mainly through radiative pumping. 

Fourteen absorption systems of molecular hydrogen 
have been identified to date in the spectra of distant quasars 
(|Noterdaeme et al. 20081 ISrianand et al. 2008j) . One of 
these belongs to the spectrum of QSO 1232+082. The pecu- 
liarity of this system is that the effective width parameter 
of the velocity distribution b for the excited rotational levels 
of hydrogen molecules is larger than the parameter b for HD 
molecules and CI atoms (Ivanchik et al., to be published). 
In this paper, we will show how the propagation of direc- 
tional radiation can broaden the velocity distribution of H2 
molecules at the levels populated by radiative pumping. In 
addition to our effect, we will consider a model that takes 
into account the presence of a photodissociation region, as 
was suggested by INoterdaeme et al. 20071 

OBSERVATIONAL DATA 



The Absorption System in the Spectrum QSO 1232+082 

The absorption system in th e spectrum of QSO 12 32+082 
has been known since 1999 (|Ge fc B cchtold 1999). It was 
measured in more detail by |Petitjean et al. 2000| with the 
8.2-m VLT telescope using theUVES echelle spectrograph. 
Subsequently, in 2001, HD mol ecules were first detected 
in this system at high redshifts ( Varsh aTovich et al. 2001|) . 
This had long been the only identification of HD molecules 
at high redshifts (the second system was identified by 
ISrianand et al. 2008|) . Analysis of this system revealed ev- 
idence for a broadening of the velocity distribution of 
molecules at excited rotational levels of molecular hydro- 
gen. This effect is important in that assigning the same 
Doppler parameter to all rotational levels leads to a sys- 
tematic error in the column density and, as a result, in 
other physical parameters of the system. 



The data on the absorption system of QSO 1232+082 
used here are presented in Fig. [2j We used the curve-of- 
growth method to analyze the data. The curve of growth 
shows the dependence of the line equivalent width on the 
column density. Since it also depends on the transition line 
parameters, the column density of an element and the dis- 
tribution parameters can be determined in certain cases 
given the set of various lines. For a Maxwellian distribution 
(with the Doppler parameter b) , the curve of growth is said 
to be a standard one. In our paper, we use a noncquilib- 
rium velocity distribution to construct the curve of growth 
and call this curve a nonstandard one. Figure shows the 
equivalent widths for various excited rotational levels of 
H2 (J = 2,3,4,5), HD, and CI. In addition, the standard 
curves of growth are plotted in the figure for b — 4.5km/s 
b = 1.77 km/s. The equivalent widths for the HD lines 
(upward triangles) fall best on the curve of growth with 
b = 1.77 km/s. The data on carbon (downward triangles) 
for this cloud have the Doppler parameter b = 1.67 km/s. 
It thus follows that the equivalent widths of the H2 lines 
should have had b equal to b = 1.82 km/s. However, the 
excited rotational levels of molecular hydrogen fall best on 
the curve of growth with b = 4.5 km/s, i.e., they have a dis- 
tribution that is two and a half times broader than that ex- 
pected from our analysis of the HD and C I data. Note that 
the Doppler parameter b for the J = 0, 1 rotational levels 
cannot be determined by the curve-of-growth method, be- 
cause the absorption lines are strongly saturated and their 
equivalent widths fall on the square-root segment where 
the curves of growth with various b merge together. Thus, 
a broadening of the velocity distribution at excited rota- 
tional levels of molecular hydrogen manifests itself in this 
system. 
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Figure 2. Standard curves of growth for the H2, HD, and CI lines corresponding to the z = 2.33771 absorption system 
of QSO 1232+082. The equivalent widths for various excited rotational H 2 levels are plotted: J = 2 (squares), J = 3 
(circles), J = 4(diamonds), and J = 5(astcrisks). The equivalent widths of the transitions from the excited rotational 
levels fall best on the standard curve of growth with b = 4.5km/s. Also shown are the HD (upward triangles) and CI 
(downward triangles) data. The standard curves of growth with b = 1.77 km/s b — 1.67 km/s are presented for HD and 
CI, respectively. 



Observations of the Broadening Effect in Other Systems 

A broadening of the effective Doppler parameter b 
with increasing rotational level J was observed pre- 
viously for our Galaxy ( Spitzer fe Cochran 1973[ 
IJenkins fc Peimbert 19971 ILacour et al. 2005|) an d 
in a high-redshift s ystem INoterdaeme et al. 2007[) . 
Spit zer fe Cochran 1973| were the first to notice this effect 
while analyzing the curve of growth for observations 
with the ultraviolet telescope of the Copernicus Orbital 
Observatory. However, since the data had large errors, it is 
hard to say whether this effect is systematic. Much later, 
the observations by IJenkins fc Peim bert 1997| on IMAPS 
revealed a broadening in molecular hydrogen toward 
£ Ori A. Several years ago, while observing excited rota- 
tional levels of molecular hydrogen in the directions of four 
young stars from the FUSE satellite, ILacour et al. 2005 
found that such an effect was present in all four direc- 
tions. The fourth and most recent observation refers to 
extragalactic molecular hydrogen. Notcrdac me et al. 2007 
revealed this effect for the absorption system at redshift 
z = 2.402 in the spectrum of the quasar HE 0027-1836. 

Note that no satisfactory explanation for the broad- 
ening effect has existed until now. Different authors use 
different models to explain it. The system observed by 
IJenkins fc Peimbert 19971 is peculiar in that the line pro- 
files exhibit three components for which the column densi- 
ties at rotational levels are low (Nj < 10 15 sm -2 and that 
there exists a systematic shift in the velocities with increas- 
ing rotational level. This enables the authors to explain the 
effect in question by a model with a shock. The observed 
broadening and the velocity shifts stem from the fact that 
the flow from the shock front gradually cools down and 



decelerates. This model runs into difficulties in explaining 
the data in which one component of molecular hydrogen is 
observed, there are no velocity shifts for various rotational 
levels, and the H 2 column densities are high. 

ILacour et al. 20051 assume that during its formation on 
dust, the hydrogen molecule acquires an additional energy 
of ~ 4.5 eV that is distributed between the vibrational, 
rotational, and translational degrees of freedom. In other 
words, the molecules are formed in excited states with an 
excess kinetic energy. However, the authors point out that 
the formation rate of molecules on dust should be an order 
of magnitude higher to explain the broadening effect. This 
is particularly problematic for the high-redshift absorption 
systems of molecular hydrogen associated with damped 
Ly-a systems; observations of the metallicity in the latter 
(|Ledoux et al. 20031 INoterdaeme et al. 2008ft show that the 
dust content in them is more than an order of magnitude 
lower than that in our Galaxy. The damped Ly-a systems 
are the systems with very high atomic hydrogen column 
densities (N(HI) > 10 20 sm~ 2 ) that are observed in quasar 
spectra and that are probably high-redshift protogalaxies. 

Another explanation for the broadening effect was also 
offered by INoterdaeme et al. 20071 They showed that the 
broadening could be explained by the presence of two 
components in the cloud. The first component is under 
the influence of strong background ultraviolet radiation, 
which increases the radiative pumping of excited rota- 
tional levels and heats up the medium through hydrogen 
dissociationthis is the so-ca lled photodissociation region 
(Hollcnb ach fc Tielens 1999)) . In the second component, the 
ultraviolet radiation is appreciably weaker; therefore, it is 
colder and the excited levels are populated sparsely. In 
that case, the excited and ground rotational levels man- 
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Figure 3. Geometric model of cloud. Directional radiation 
with an intensity z) falls on the cloud. 

ifest themselves mainly in the hot and cold components, 
respectively, which can provide the observed effect. As will 
be shown below, when the radiative transfer in the Lyman 
and Werner bands is considered, these regions naturally 
emerge in the same molecular cloud. 

Note that ILacour et al. 20051 INoterdaeme et al. 20071 
assumed that radiative pumping could not produce a broad- 
ening. The authors placed emphasis on the fact that the 
molecule did not change its velocity during radiative pump- 
ing transitions and sought for other possible explanations. 
This assertion is beyond doubt. However, we will show in 
the next section that when the transfer of directional ra- 
diation is considered in detail, a broadening at the levels 
being populated can arise during radiative pumping if the 
population from the line wings is taken into account. 

MAIN POINTS OF THE MODEL 

A broadening of the velocity distribution at the levels 
populated by radiative pumping can be obtained by consid- 
ering the radiative transfer. The broadening effect will be 
at a maximum during the propagation of directional radia- 
tion. Therefore, let us consider a molecular hydrogen cloud 
on which directional radiation with an intensity z) falls 
(see Fig. [3]) . We will be interested in the velocity distribu- 
tion of the molecules populated by radiative pumping as 
a function of the depth of radiation penetration into the 
cloud. 

To understand the effect qualitatively, let us consider a 
three- level system (see Fig. 0}. Naturally, this system will 
fully reflect the main physical properties of the hydrogen 
molecule of interest to us and it can be used as an example 
to obtain the broadening effect analytically. Thus, we have 
the ground level (0) that acts as a reservoir. A metastable 
level I is radiatively pumped from the ground level through 
the systems photoexcitations to level k followed by its re- 
laxation to level I. The direct transitions from the ground 
level to level I will be assumed to be forbidden. 

Based on the characteristic time scales shown in Fig. [4] 
(which correspond to those for the hydrogen molecule) , we 
can assume that the system passes from level k to level I 
almost instantly. Since level / is metastable, we can observe 
the molecules at this level. To obtain the velocity distribu- 
tion of the molecules at level I, let us consider the balance 
equation for the molecules at the level 

ni(v Zl z)Ai = n (v z ,z)P ok (v z ,z) (1) 

The left-hand and right-hand sides of the equation describe, 
respectively, the spontaneous downward transitions and the 
population of level I from the ground level through radiative 
pumping via level k. Pofe is the capability of the particles to 
pass to level k under radiation or the fraction of the parti- 
cles that pass to level k per unit time at depth z and that 
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Figure 4. Three-level system: is the ground level from 
which level / is radiatively pumped through the resonance 
level k. The characteristic time scales of the processes con- 
sidered in this system corresponding to a hydrogen molecule 
in the interstellar medium are shown near the arrows indi- 
cating these processes. 

have a velocity v z ; in other words, Pofc is the photoexcita- 
tion rate: 

+ oo 

/I (is z\ V 
^'J K 0k (v + v^j dv. (2) 

— oo 

It is important to note that this quantity is a function of 

(Vz,z). 

We can neglect the multiple scatterings while consider- 
ing the radiative transfer in the cloud, since the resonance 
photons that accomplish radiative pumping fragment as the 
excited states decay due to the presence of a ro- vibrational 
cascade in the H2 molecule. Therefore, we can use a radia- 
tive transfer equation in the form 

I(y,z) ^ I(iy,0)e- Tok ^ z \ (3) 

where 

z +00 

T 0k {v,z) = J dz' J dv z n {v z ,z')n ok (v + . (4) 
-00 

The quantity «ofc in Eqs. and (UJ) is the photoexcitation 
cross section for the — > k transition. 

The Broadening Effect 

An expression for the velocity distribution of the molecules 
at level I can be obtained by rewriting Eq. |T]) using Eq. 

n Q (v z ,z) f I(u,z) f v z \ 
n l (v z ,z) = — / — n ok [v + v— j dv. (5) 

—00 

The distribution at level I is broadened during the forma- 
tion of an absorption line in the — > k transition followed 
by the population of level I (Fig. [4]). Figure [5] presents the 
calculated number and column densities and line profiles 
for typical molecular hydrogen line parameters (we chose 
the L3-0R(0) line). In what follows, we assume that the 
distribution at the ground level is Maxwellian. 

The upper panels in Fig.[5](al, bl, cl) present the line 
profile for a given optical depth ro at the line center (indi- 
cated above the panels). The middle panels (a2, b2, c2) 
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Figure 5. Spectral resonance line profiles and velocity distributions of molecules. The calculation is presented for three 
optical depths at the line center: (a) r = 0.2, (b) r = 20, and (c) r = 1000. The upper panels show the profile of the 
line through which level / is radiatively populated. The middle panels show the normalized molecular number density 
distribution produced by radiative pumping. The lower panels show the normalized molecular column density distribution 
integrated along the line of sight to the depth corresponding to the line profile on the upper panels. 



show the normalized distribution of the number density 
populated by radiative pumping to level I at depth z where 
the line profile was formed. The presented number density 
distributions were normalized to the area under the distri- 
bution for To = 0.2 (panel (a2)). Therefore, the area under 
the distributions for r = 20 (panel (b2)) and r = 1000 
(panel (c2)) shows the extent to which the total photoexci- 
tation rate decreases compared to the depth where r = 0.2. 
The lower panels (a3, b3, c3) display the normalized dis- 
tribution of the column density, i.e., integrated over the 
line of sight to depth z. The column density distributions 
were normalized to the area of the distribution under its 
curve for To = 1000 (panel (c3)). The velocity in km/s is 
along the x axes (for the line profile, this corresponds to 
the Doppler wavelength shift). We see that as long as the 
line profile is unsaturated (t = 0.2, (al)), the populated 
(a2) and column (a3) densities correspond to a Maxwellian 
distribution at the ground level, i.e., no broadening occurs. 
When the lines are saturated (to = 20, (bl)), the distri- 
bution populated to level I becomes a highly nonequilib- 
rium one. As long as the absorption is a Doppler one, the 
molecule absorbs a photon that closely corresponds to its 
Doppler velocity shift. In other words, the photoexcitation 
cross section K a k in Eq. §5§ can be replaced with a delta 
function; therefore, the dependence of the photoexcitation 
cross section P k on the velocity of the molecule will closely 
correspond to the profile of the — * k transition line. When 
the line is saturated, the photons are deficient at the line 
center and the population takes place mainly at the line 



edges, where the photons are plenty. Therefore, the num- 
ber density distribution at level / becomes a highly nonequi- 
librium one (panel (b2)). The column density distribution 
(panel (b3)) is an observable quantity in the spectra; it ex- 
hibits a broadening with the formation of a plateau. 

As the line is saturated further (r = 1000, (cl)), i.e., 
as the radiation penetrates deeper into the medium, the 
Lorentz wings will manifest themselves in the line profile. 
This means that level I will also be populated in the Lorentz 
wings. In this case, the transition wavelength does not cor- 
responds to the Doppler velocity shift and the replacement 
of «ofc with a delta function is invalid. This population is 
peculiar in that the distribution populated by the Lorentz 
wings will correspond to the distribution at the ground 
level, i.e., it will be Gaussian. We see from panel (c2) that 
even a small fraction at a given optical depth corresponds 
to a noncquilibrium distribution. Thus, the population in 
the Lorentz wings can greatly reduce the broadening of the 
column density distribution (panel (c3)). 

The nonequilibrium distribution populated to level I will 
manifest itself in the absorption lines from this level, while 
the lines will not be described by a Voigt profile. When a 
Voigt profile is fitted into the line, the effective width pa- 
rameter of the Maxwellian distribution b will naturally be 
larger than that for the ground level. Therefore, when we 
analyze the absorption lines for level I and use the stan- 
dard theory of radiative transfer, the parameter b to be 
determined will effectively increase. 
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Thermalization and Turbulence 

The nonequilibrium distribution populated to level I can be 
thermalized with time or remain broadened in the presence 
of turbulent motions. It all depends on the relationships 
between the rates of characteristic processes. To describe 
the turbulent motions, we will use the assumption about 
microturbulence under which the distribution in the cloud 
in turbulent velocity w tur b is described by a Gaussian func- 
tion with a parameter bturb- The ground level may be con- 
sidered to be thermalized, i.e., the velocity distribution of 
the molecules related to thermal motion, w t h, is described 
by a Maxwellian distribution with a parameter bth- The to- 
tal distribution obtained by a convolution of the turbulent 
and thermal dist ributions w ill then be Maxwellian with the 
parameter b = \/ b 2 urb + b\ h Therefore, the relationship be- 
tween the following three rates of the processes should be 
considered: the thermalization rate, the turbulization rate, 
and the photoexcitation rate from level The turbuliza- 
tion rate may be considered to be much lower than the 
other two rates. If the thermalization rate is higher than 
the photoexcitation rate from level I, then the distribution 
that manifests itself in the absorption line from level / will 
be thermalized. 

The distribution thermalization effect for static turbu- 
lence can be written as 



hi (v z , z) = / dv' z ni (v' z ,z)w + 



b 2 

Z b 2 



(6) 



(for the derivation, see the Appendix). 

When bth — * w+(v) — > 6(v) and, hence, hi(v,z) = 
ni(v,z), i.e., the nonequilibrium distribution is retained. 
When bturb —> 0, w + does not depend on v z and the ve- 
locity distribution is completely thermalized, acquiring a 
Gaussian shape with b = b t h- 



Optical Depths of Broadened Lines 

As a result, we find that, given thermalization ([B]), the dis- 
tribution at level I will be 



©• Thus, we obtain 

L +oc 
Ti m (v,L) = J dz 

-oo 
L +oc 



dv z K tm (u + v^j hi(v z ,z) 



= J dz J dv z K lm [v + v—j x 

-oo 
+ oo 

x J dv'Mv' z ,z)w + L z -v' z b ^\ 



(8) 



It is important to note that in our approximations, the 
z dependence is present only in ni(v z , z). Using Eq. (O, we 
can represent the expression for the optical depth as 



Tlm{v) 



(9) 



where Wi{v z ) is the nonequilibrium distribution of the par- 
ticles thrown to level I integrated along the line of sight: 



Wi(v z ) 



Ni(v z ) 



J dz ni(v z ,z) 



+oo L +oo 

J dv z Ni(v z ) J dz J dv z ni(v z ,z) 

— oo — oo 



(10) 



The quantity T^ m in Eq. © is the optical depth in the 
standard theory of radiative transfer, which can be written 
via the Voigt function as 



v' h 2 
c b z 



- e 2 AntNifi 

: 

m e c 7 



where 

air, 

x - 



■H(ai m ,x), 
(11) 



7c 



47W /m b th \/T+ {hurb /b) 2 

c , v 'z b turb v Vn 



bthV 1 + ( b turb/by 



b 2 



We can write expressions for Ni(v' z ) the column density 
distribution of the particles that are populated to level I 
and that have no time to be thermalized. Using Eq. ([5]), we 
will find that 



+ OC 



ni{v z ,z) = 1 ^f!- / dv z P ok (v' z ,z)x 



A, 



x n (v' z , z)w+ I v 



b 2 

I u turb 
Z b 2 



(7) 



Ni(v z ) = J dzni(v Zl z) = -i- J dzn (v z ,z)P k(v z ,z). 
o o 

(12) 

For a uniform and equilibrium distribution at the 
ground level and for the radiative transfer equation written 
for the i — > k lines in form ([3]), we can easily obtain 



As we showed above, the particle velocity distribution 
at level I §5§ can manifest itself in the absorption line asso- 
ciated with the transition from this level to some resonance 
level m. We will assume that the absorption line profile is 
defined by the optical depth in accordance with Eq. 
The optical depth r; m for the I m transition will have 
the form (j4]) with the cross section K; m corresponding to 
the transition and the number density hi(v z ,z) given in 



Ni{v z ) 



No 
A, 



w(v z ) 



dv x 
hv 



(1 - e-^o^i)) 



Kofc (y + , (13) 



where Iq{v) = I(v, z = 0). 

Passing to the limit in Eq. ([9]) allows physically mean- 
ingful expressions to be derived for the extreme cases: 
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1. The absence of turbulence {b tur b <C bth)- 
In this case, b t h ~ b. Then, 



+00 



dv'Mv'Jr+iv) 



'lm 



M, (14) 



where Tvt, is the optical depth in the standard theory 
of radiative transfer with a thermal broadening b = b t h- 
Thus, in this case, the nonequilibrium distribution ra- 
diatively populated upward is completely thermalized to 
become an equilibrium one and gives a standard shape 
for the line profile. 
2. Dominant turbulence {bturb 3> bth)- 



In this case, b 



turb 



b and then I i> 



Vlr, 



62 



will correspond to the cross section m m (y + v^r) to 
constant coefficients and the broadening effect related 
to radiative pumping will manifest itself in full: 



r lm (v, L) = / dv z Ki m (u + v— j Ni(v z ). (15) 



Nonstandard Curves of Growth 

Using the above formulas, we can demonstrate how the 
broadening effect will appear in curve of growth represen- 
tation. In Fig. [SI the dashed lines indicate the standard 
curves of growth for distribution parameters b from 1.8 to 
3.8 km/s at steps of 0.5 km/s. In the same figure, the circles 
indicate the curves of growth for the broadening effect that 
we calculated based on the three-level model. The curves of 
growth correspond to the transitions from level /, i.e., to the 
particle velocity distribution at level /. Since, as was shown 
above, this distribution can differ from the Maxwellian one, 
the curves of growth will also differ from the standard ones. 
The figure shows two limits: with (open circles, see ©) and 
without (filled circles) thermalization. 

The shape of the nonstandard curve of growth depends 
on certain parameters of the calculation. First, these in- 
clude the parameter b, the width of the Maxwellian dis- 
tribution at the ground level from which the population 
occurs. For the calculation in question, btot is b to t km/s, 
which corresponds to the lower standard curve of growth 
in Fig. [5] and agrees with the HD and CI data. For a ther- 
malized distribution at level I (open circles), the relation- 
ship between bth bturb should be specified. The broaden- 
ing effect will depend on this relationship. Irrespective of 
the relationship between the characteristic time scales, the 
case without thermalization (filled circles) can take place 
for bturb b t h- In this case, the curve of growth passes 
higher, i.e., it corresponds to larger effective b for the stan- 
dard curves of growth. Naturally, the limit b tur b <C bth 
corresponds to complete thermalization, with the nonstan- 
dard curve of growth coinciding with the standard one for 
b = 1.8 km/s. Second, the shape of the curve of growth de- 
pends on the column density at the ground level. This is 
easy to explain, since the shape of the distribution (Fig. [5l 
panels (c)) depends on the column density at the ground 
level in the cloud via the line profiles (Fig. [5J panels (a)). 
In the case under consideration, the column density was 
chosen to be No = 10 19 cm -2 . 



It should be noted that the broadening effect manifests 
itself only in a certain range of column densities at the 
ground level. The distribution is not yet broadened at low 
column densities {N to t < 10 15 cm -2 for typical line parame- 
ters of the Lyman or Werner molecular hydrogen band). At 
high column densities {Ntot > 10 20 cm~ 2 ), the broadening 
effect is suppressed by the Lorentz wings (see Fig. [5]). 

Thus, the inferred effective broadening takes place in 
the case of directional radiation in the cloud and using the 
standard curve of growth can lead to a systematic error in 
the column densities and distribution parameters. 



A SCHEME FOR CALCULATING A MOLECULAR 
HYDROGEN CLOUD 

Let us apply the broadening effect that we obtained in the 
previous section to molecular hydrogen. This effect arises 
when considering the transfer of directional radiation in a 
medium and is related to the peculiarity of the radiation 
interaction with the energy level structure. Therefore, we 
will consider a cloud composed only of molecular hydro- 
gen. Since the radiation directionality is important for the 
manifestation of the broadening effect, we will consider a 
one-dimensional cloud model, i.e., one direction in which 
the radiation propagates coincident with our line of sight. 
In this approximation, we neglect the background radia- 
tion compared to the directional one, which is quite valid 
when a bright (in the ultraviolet) star (or a quasar) is close 
to the cloud. Naturally, in this case, the star lies behind 
the cloud and the line of sight coincides with the direc- 
tion of the star. It is important to note that in this cloud 
model, we do not consider such processes and the photodis- 
sociation of molecular hydrogen and its formation on dust 
(jHollenbach et al. 19711 I Jura 1975)) . A remark about the 
influence of dust will be made below. 



Multilevel Model 

When we pass from the three-level system to the real mul- 
tilevel model of the hydrogen molecule, a number of pecu- 
liarities arise that should be taken into account. 

The first important peculiarity is that there exist two 
different states of molecular hydrogen - ortho- and parahy- 
drogen. In general, the direct radiative transitions between 
these states are forbidden (t~5x 10 12 yr). A much faster 
mixing of the ortho- and para-states is possible through 
collisions with H + and (Abgr aTTet al. 19~92| and dur- 
ing the dissociation of hydrogen followed by its formation. 
However, under our model assumptions, the ortho- and 
parahydrogen systems may be considered separately. 

The second peculiarity consists in an allowance for 
the presence of many lines in the Lyman and Werner 
bands through which radiative pumping proceeds. This 
means that the sum over all the possible transitions of the 
Lyman and Werner bands will now appear in Eq. @ for 
the photoexcitation rate. To calculate the latter, we used 
the wavelengths and oscillator strengths for the Lyman 
and Werner band transitions from |Abgrall et al. 19 93a, 
Abg rall et al. 1993b respectively. The weight factor that 
describes the ro-vibrational cascade in the ground elec- 
tronic state ( |Dalgarno Ik Black 1974[ ) should also be added. 
The ro-vibrational cascade is peculiar in that it allows high- 
lying rotational levels to be populated from low-lying ones 
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Figure 6. Nonstandard curves of growth. The circles indicate the curves of growth for the velocity distribution of the 
molecules broadened by radiative pumping: the filled and open circles indicate the results of our calculation, respectively, 
without and with thermalization of the velocity distribution at the level. The following distribution parameters were 
chosen: b t h = 0.8km/s, b tur b — 1.6km/s. The dashed lines represent the standard curves of growth with parameters 
b from 1.8 to 3.8 km/s at steps of 0.5 km/s. In the absence of a broadening effect, the nonstandard curves of growth 
marked by the circles would coincide with the lower standard curve of growth with b — 1.8 km/s. 



(jJura 1975|) . We will assume that the ro-vibrational cascade 
occurs almost instantly, i.e., on time scales much shorter 
than those of the processes considered in our model. 

The next remark is that the first several rotational levels 
of the hydrogen molecule should be considered jointly, i.e., 
we should depart from the consideration of the ground level 
as a reservoir. For example, the levels only with even rota- 
tional numbers (J — 0, 2, 4, etc.) are present in the parahy- 
drogen system in the ground electronic state. Assuming 
that almost all of the molecules are at the lower vibrational 
level, we should take into account the radiative pump- 
ing of the J = 2 level from the J = level and vice 
versa. We can restrict our analysis to the first several ro- 
tational levels, since the time scales of the spontaneous 
J — > J — 2 transitions in the lower vibrational state de- 
crease with increasing J; this means that the upper rota- 
tional levels can be included in the ro-vibrational cascade. 
The data for the ro-vibrational cascade that we used from 
|Dalgarno fc Black 1974| are restricted to the first eleven ro- 
tational levels, which are enough under typical conditions 
of the interstellar medium. 

Yet another peculiarity is that although we disregarded 
the presence of atomic hydrogen, it can be assumed that 
there is a dip in the spectrum for A < 912 A related to the 
absorption of atomic hydrogen in the ionization continuum 
for the radiation incident on the cloud. This allows us to 
restrict our analysis of the transitions in the Lyman and 
Werner bands to the 20th and 6th vibrational levels, re- 
spectively. As a result, we calculated the radiative transfer 
for approximately 600 resonance lines of the Lyman and 
Werner H2 bands. 



Balance Equations 

Let us write the balance equations for the system of molec- 
ular hydrogen levels under our assumptions. In a stationary 
situation, 

dn^Vz, z) _ 

at 

and it is convenient to write the system of equations to 
determine the particle distributions at the levels in matrix 
form, 

G(v z , z)n(v z , z) = 0, (16) 

where n(v z , z) is the column of the particle velocity distri- 
butions for various rotational levels; G is the matrix that 
describes the transitions between the levels, 

G(v z ,z) = B(v z ,z)+A + C, (17) 

where the matrix A is related to spontaneous quadrupole 
transitions between rotational states, the matrix C is re- 
lated to collisional transitions, and the matrix B is related 
to radiative pumping. The specific form of the matrix can 
be easily obtained by considering the balance between the 
particles going from and coming to some level i. For ex- 
ample, the components of the matrix B can be represented 
as 

Bij =Uij -^ ^n ife , (18) 

k 

where Sij is the Kronecker symbol and Ily is the photoexci- 
tation rate for the transition from the zth rotational level to 
the jth one and, what is important, is a function of (v z , z) 
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and depends on rii(v z , z). This is related to radiative trans- 
fer, which manifests itself in the saturation of the radiative 
pumping line, and it is this fact that provides the required 
broadening. The photoexcitation rates can be written as 



(19) 



where k is the set (B(C), v' , J') (J' are the permitted rota- 
tional levels). Accordingly, Pik is the photoexcitation rate 
for the i — > k transition: 



+00 



Pik(v z ,z) 



i (y, f) 

hv 



dv. 



and Kkj is the probability to relax k =(B(C), u', J 1 ) — + (X, 
v = 0, J = i); for convenience, it can be separated into 
the relaxation from excited electronic levels and the ro- 
vibrational cascade, 



J)i ■ 



(20) 



where G describes the relaxation from upper levels and 
(3 is responsible for the cascade in the lower electronic 
state. We used the coefficients of the rovibrational cas- 
cade from Dalgarno & Black 1974 The collisional transi- 
tion rates (matrix C) were taken from[Abgrall ct al. 1992 
System (ITTJl) is underdetermined, since one of the equa- 
tions is not linearly independent. Therefore, we should ex- 
clude it and expand the system by an additional natural 
conditionthe condition on the total number density. Since, 
as we noted above, the ortho- and parahydrogen systems 
may be considered independently, the conditions on their 
total number density can be written separately: 

no(v z ,z) + n 2 {v z ,z) + s + n w {v z ,z) = n para (v z ,z), (21) 

for parahydrogen and a similar condition for orthohydro- 
gen (with the total number density n ort ho)- n para (v z , z) 
and n ort ho{v z , z) are model-dependent quantities. We per- 
formed the calculation by assuming that these distributions 
were Maxwellian and that the integral of n para (v z , z) + 
n or tho{v z , z) over all velocities was equal to the total num- 
ber density of molecular hydrogen. Assuming that these 
systems are in thermodynamic equilibrium, the ratio of the 
para- and orthohydrogen number densities is determined 
by the temperature in the cloud (typically, T ss 100° K). 

System (fT!)]) should be solved numerically. We consider 
a homogeneous medium with typical molecular hydrogen 
number densities, nn 2 = 10-7- 1000 ~ 3 . N tot (z), the column 
density of molecular hydrogen to depth z, can serve as a 
measure of the depth of penetration into the cloud. 



RESULTS 



Main Results 

The results of our calculations based on the scheme pre- 
sented in the previous section are shown in Fig. [7J for 
various combinations of N to t (the total molecular hydro- 
gen column density in the cloud) and Iq (the intensity of 



the incident directional radiation). For two column densi- 
ties Ntot, 10 18 cm -2 10 19 cm -2 , we chose three different 
intensities of the incident radiation: 10 -16 , 10 -17 10 -18 
erg/cm 2 /s/Hz. Thus, there are a total of six different cal- 
culations (Fig. [7]). The background isotropic radiation was 
assumed to be negligible compared to the directional one 
(this is valid when the directional radiation is determined 
by a nearby star or a quasar). For all calculations, the dis- 
tribution parameters were chosen from observational data 
on HD and CI in the QSO 1232+082 absorption sys- 
tem: bth — 0.81 km/s for the thermal distribution and 
bturb — 1.6 km/s for the turbulent one; as a result, the 
Doppler parameter b for the distributions n para (v z ) and 

n or to{v z ) (see flUD is b = y/bj h + b'f urb = 1.8 km/s. The 
total molecular hydrogen number density was assumed to 
be uniform and the calculation was performed for several 



values of Uq 



L para 



10, 100, 1000 enr 



As the result of our calculation, Fig. [7J presents the de- 
pendences of the line profiles for the Lyman and Werner 
molecular hydrogen bands on the depth of radiation pen- 
etration into the cloud. The line profiles reflect the parti- 
cle velocity distribution on the line of sight at these levels 
and are calculated consistently. To present our results, we 
chose several lines for which the equivalent width was cal- 
culated and plotted on the curve of growth. For compari- 
son with the observational data for the absorption system 
of QSO 1232+082 , we used the transitions from the rota- 
tional J = 2(squares), J = 3(circles), and J = 4 (asterisks) 
levels. Figure [7J present two extreme cases: the absence of 
thermalization (filled symbols) and complete thermaliza- 
tion at all levels (open symbols) for the chosen parameters 

bth bturb- 

As was shown above, one might expect a significant 
broadening of the particle velocity distribution (by almost 
a factor of 2 in terms of the effective Doppler parameter) 
in the three-level system. However, the broadening of the 
distribution in the cloud model decreases appreciably. The 
main reason is that there are many radiative population 
transitions. For the chosen rotational level, there are about 
60 such transitions in the Lyman or Werner bands. Each of 
these transitions has its own values of A/ the product of the 
transition wavelength by the oscillator strength. For each of 
them, the broadening effect will begin (with line saturation) 
and be suppressed (when the Lorentz wings appear) at dif- 
ferent depths in the cloud. This is easy to understand, since 
the optical depth includes the combination NXf. Therefore, 
the incompletely broadened distributions from some transi- 
tions and the distributions with the influence of the Lorentz 
wings (reducing the broadening) from other transitions will 
be admixed to the maximum broadening of the distribution 
from some transition at some depth. 

In addition to the aforesaid, the situation will be com- 
plicated at a high intensity of the incident radiation. In 
this case, the excited rotational levels will be populated 
more heavily than the ground levels until the lines of var- 
ious transitions become saturated and the photoexcitation 
rate decreases. This will be true for the cloud shell. We will 
consider the peculiarity that arises in this case below. 

When the curve of growth is analyzed, it is necessary 
that the lines from the levels with a broadening fall on the 
logarithmic segment of the curve of growth. Otherwise, if 
the broadening is present, we will not see it on the linear 
and square-root segments due to the peculiarity of the curve 
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Figure 7. Results of our calculation of the molecular hydrogen cloud model for various intensities of the incident radi- 
ation I and total molecular hydrogen column densities N to t- The squares, circles, and asterisks indicate the equivalent 
widths of the lines from the J = 2, 3, 4 levels, respectively. The filled and open symbols represent the results of our calcu- 
lation, respectively, without and with thermalization of the velocity distribution at the level. The following distribution 
parameters were chosen: b t h = 0.8 km/s, b tur b — 1.6 km/s. The dashed lines represent the standard curves of growth with 
parameters b from 1.8 to 3.8 km/s at steps of 0.5 km/s. In the absence of a broadening effect, all symbols would fall on 
the lower standard curve of growth with b = 1.8 km/s. 
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of growth. Therefore, when the system of molecular hydro- 
gen is considered, the dependence in the manifestation of 
the broadening effect on the intensity of the radiation inci- 
dent on the cloud is obvious. The excited rotational levels 
will be populated more heavily with increasing intensity 
of the incident radiation. This is clearly seen if we com- 
pare panel (a) in Fig. [7] with panels (c) and (e) (and panel 
(b) with panels (d) and (f)). In this case, the points corre- 
sponding to the transitions from different rotational levels 
become closer to one another (cf. the squares, circles, and 
asterisks). 

The broadening effect in the chosen six combinations 
of parameters reaches its maximum at N to t = 10 18 cm~ 2 , 
la = 1CP 17 erg/cm 2 /s/Hz (panel (c) in Fig. [7]). This can 
be explained as follows: the total hydrogen column density 
N to t should not be high (> 10 19 ~ 2 ); otherwise, the Lorcntz 
wings begin to suppress the effect (see Figs.[7p,[7p!, and[?F). 
Obviously, the radiation must be intense enough for the 
excited rotational levels to fall into the region where the 
logarithmic part of the curve of growth begins, where the 
effect is most pronounced. If the intensity is higher, then 
the excited levels are populated more heavily (Fig. [7fe) and 
approach the square-root part of the curve of growth; if the 
intensity is lower (Fig. [7K), then the levels fall on the linear 
part of the curve of growth. On these segments of the curve 
of growth, the broadening effect does not manifest itself, 
because it is insensitive to b. 

However, the transfer of directional radiation that we 
considered is not enough to explain the broadening of 
the distribution for excited rotational levels of molec- 
ular hydrogen that arises in the absorption system of 
QSO 1232+082. Therefore, additional possibilities for the 
explanation should be sought. Nevertheless, this effect can- 
not be neglected, provided that the radiation directionality 
in the cloud is dominant. 



The Influence of Dust 

One of the ways to explain the broadening of the velocity 
distribution is the formation of molecular hydrogen on dust 
(|Lacour et al. 2005]) . However, the dust-related broadening 
effect during the propagation of directional radiation can 
be only a small correction. When a stationary situation is 
considered, the amount of formed hydrogen is assumed to 
be exactly equal to the amount of dissociated one, which 
accounts for 13% of the molecules excited in the Lyman or 
Werner bands. The remaining 87% are involved in radiative 
pumping, which, as was shown above, leads to a broaden- 
ing of the distribution. This means that the two effects will 
be added with weights of 0.87 and 0.13 for the broadenings 
caused by radiative pumping and dust formation, respec- 
tively, relative to the distribution functions. 

The presence of dust in the cloud can lead to a slightly 
different effect. The point is that dust reduces the intensity 
in the continuum. In fact, the population by the Lorentz 
wings also takes place in the continuum. Therefore, in the 
presence of dust in the cloud, it is capable of eating away 
the ultraviolet continuum in such a way that starting from 
some hydrogen column density Nd (depending on the dust 
number density in the cloud), it essentially will reduce the 
population in the Lorentz wings to zero, i.e., in an ideal 
case, it will only freeze the broadening. We considered this 



additional effect, but it will not allow us to obtain the addi- 
tional broadening needed to explain the observational data. 

The Model of a Hot Shell 

The assertion of [Notcrdae me et al. 20071 that the velocity 
distribution at excited rotational levels can be broadened in 
the presence of a photodissociation region seems the most 
realistic explanation. This explanation naturally also arose 
in our paper when we considered the results of our model 
calculation presented in Fig. [5) 

In Fig. [SJ the populations of various rotational lev- 
els are plotted against the depth of radiation penetration 
into the cloud. We see from the figure that the cloud can 
be arbitrarily separated into two regions: the shell and 
the inner part. In the shell, the lines of the Lyman and 
Werner bands are not saturated and the radiative pump- 
ing rate is high at a sufficiently high intensity (in our 
case, Iq = 10~ 17 ers/cm /s/Hz). Consequently, the transi- 
tion rates between the rotational levels of molecular hydro- 
gen in the lower vibrational (ground) electronic state are 
also high. In that case, the excited rotational levels will 
be populated significantly (the J = 2, 3 and 4 levels are 
populated more heavily than the ground ones). Since the 
absorption lines become saturated as the radiation pene- 
trates into the cloud (starting from Ntot ~ 10 15 cm -2 ), the 
photoexcitation rates and, hence, the population of the ex- 
cited rotational levels decrease. As a result, only the ground 
levels turn out to be populated inside the cloud (in our case, 
Not ^ 10 16 cm -2 ) We also see that the inner part of a cloud 
with Not > 10 18 cm -2 (measured in the molecular hydro- 
gen column density) will be several orders of magnitude 
larger than its shell. Therefore, the ground levels manifest 
themselves and allow us to determine the physical condi- 
tions only in the inner part of the cloud. As regards the 
excited levels, the column density of the excited rotational 
levels in the shell and the inner part of the cloud may turn 
out to be comparable in view of the radiative population 
pattern. This means that the excited rotational levels will 
reflect the physical conditions in the shell. 

Since the photoexcitation rates in the shell are 
higher, the photodissociation rates are also higher. 
Photodissociation can heat the medium, because an energy 
of 4, 48 eV is released in one act of this process. Therefore, 
in the regions with active photodissociation, the medium 
can heat up and, as a result, have a broader velocity dis- 
tribution of the molecules. These are called photodissocia- 
tion regions (PDRs) (see, e.g.. iHollcnba ch fc Tielens 1999|, 
ILe Petit et al. 20061) . 

We will take into account the PDR by specifying a 
nonuniformity of the Doppler parameter b in the cloud. 
The medium in the shell will have a higher value of b than 
that inside the cloud. In general, b can increase in the shell 
not only through an increase in b t h, i-e., the heating of the 
medium, but also through an increase in the turbulent ve- 
locity component b tur b- To specify the gradient in 6, we will 
take into account the fact that the photodissociation rate 
(and, hence, the heating of the medium) is related to the 
photoexcitation rate. The dependence of b on the depth of 
penetration into the cloud that we chose is shown in Fig. [51 
We chose the characteristic depth at which the shellinner 
part transition occurs to be N to t = 3 x 10 16 cm~ 2 . 
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Figure 8. Populations of various rotational levels versus depth of penetration into the cloud. The cloud is arbitrarily 
separated by the vertical solid line into the shell (photo-dissociation region) and the inner part. The dotted vertical lines 
mark the depths at which the J — 1 and J — 3 (orthohydrogen) populations and the J = and J — 2 (parahydrogen) 
populations become equal. 
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Figure 9. Results of our calculation of the molecular cloud model in the presence of a hot shell introduced through a 
nonuniformity of b over the cloud (shown in the insert). The filled, open, and crossed symbols indicate the results of 
our calculation, respectively, without and with the shell and with a nonuniformity of b but without any nonequilibrium 
radiative population. The dashed lines indicate the standard curves of growth with parameters b from 1.8 to 3.8 km/s at 
steps of 0.5 km/s. 



Figure [5] presents the results of our calculation of the 
model for radiative transfer in a molecular hydrogen cloud 
with a PDR specified by a nonuniformity of the parameter 
b. The open symbols indicate the equivalent widths of the 
lines from various rotational H2 levels obtained without in- 
cluding the PDR but in the presence of a nonequilibrium 
radiative population (this calculation was presented above, 



in Fig. [71 panel (d)). The filled symbols represent the results 
of our calculation when adding a b nonuniformity shown in 
the insert in Fig. [9[ We see that the line equivalent widths 
for the transitions from excited rotational levels in the pres- 
ence of a & nonuniformity lie high, near the curves of growth 
with b ~ 4.0 km/s. This result shows that the presence of 
a PDR can be used to explain the observational data for 
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Figure 10. Possible geometric configurations of PDRs relative to the line of sight: (a) the line of sight is directed toward 
the source (a young star, a quasar) and, therefore, a PDR can be observed; in this case, the broadening effect should be 
taken into account; (b) the line of sight is perpendicular to the direction of the star and the manifestation of a PDR in 
observations is unlikely; in this case, the broadening effect is suppressed. 



the absorption system of QSO 1232+082. Because of the 
simplified scheme of our calculation, we did not attempt to 
determine the exact parameters of the nonuniformity but 
only show the importance of taking into account the PDR. 

We also see from Fig. |H] that when applying an analy- 
sis of the standard curve of growth to the computed points 
in the presence of a PDR (filled symbols), to achieve cor- 
respondence to some standard curve of growth, it will be 
necessary to displace the points greatly leftward. This will 
lead to an even larger increase in the effective parameter 
b and larger errors in the populations of excited rotational 
levels. Therefore, it is important to note that assigning the 
same Doppler parameter to all rotational levels of molec- 
ular hydrogen when analyzing the absorption systems can 
lead to serious systematic errors in the column density and 
other physical parameters. 

The crossed symbols indicate the results of our calcula- 
tion in the absence of any noncquilibrium radiative popu- 
lation, i.e., the velocity distributions at all rotational levels 
were assumed to correspond to a Maxwellian distribution 
with a nonuniform b over the cloud. We see that the ef- 
fect should be taken into account for a detailed analysis of 
PDRs. If we associate the presence of a PDR with the loca- 
tion of the cloud near a young bright star, then this effect 
should be taken into account, because, in this case, the ra- 
diation can be assumed to be directional, i.e., to propagate 
near the line of sight. The point is that a PDR close to a 
star must be located in the cloud on the side of the star 
(see Fig. fTUk). In the opposite case, where the young star- 
cloud line is perpendicular to the line of sight, the effect 
of a nonequilibrium radiative population will be unobserv- 
able on the cloud. However, in this geometry, the PDR is 
also difficult to observe, because the molecular hydrogen 
column density is low (No < 10 15 cm -2 ) in this region (see 
Fig. [TUb). In this case, if the line of sight passes through 
the PDR, then it will not pass through the cloud interior, 
where the bulk of the molecular hydrogen is concentrated, 
and, conversely, if the line of sight passes through the inte- 
rior, it will not touch the PDR. 

CONCLUSIONS 

Observations of molecular hydrogen clouds reveal a 
broadening of the velocity distribution at excited rota- 
tional levels. This effect is present in the observational data 
for Galactic ([Jenkins fc Peimbert 1997|,ILacour et al. 2005]) 



and high-redshift (jNoterdaeme et al. 2007[) clouds. We sug- 
gested a mechanism that could be responsible for the ef- 
fective broadening of the distribution at excited rotational 
levels. The condition for the manifestation of this effect is 
radiation directionality. Using a three-level model system, 
we showed that the radiative pumping of an excited level 
for directional radiation could lead to a broadening of the 
particle velocity distribution at this level. The broadening 
of the distribution was considered in two limits: with ther- 
malization of the velocity distribution (the turbulent com- 
ponent is retained) and without its thermalization. In the 
absence of thermalization, our mechanism in the three-level 
system considered (at transition line parameters typical of 
molecular hydrogen) can lead to an almost twofold effec- 
tive broadening of the distribution in an analysis using the 
curve-of-growth method. The suppression of the effect dur- 
ing thermalization depends on the ratio of bturb and bth, 
the parameters of the Maxwellian distribution attributed 
to turbulence and thermal motions 

We applied the inferred effect to the system of molecular 
hydrogen energy levels and calculated a simple molecular 
hydrogen cloud model. We considered the balance at the 
excited rotational levels of the lower vibrational ground 
electronic state by taking into account the set of lines of 
the Lyman and Werner bands through which the radiative 
pumping proceeds. We also took into account the presence 
of a ro-vibrational cascade in the ground electronic state. 
The large number of lines slightly reduces the effect, be- 
cause its manifestation in the radiative level population 
is heterogeneous. However, the broadening of the veloc- 
ity distribution of molecules during radiative pumping is 
not enough to explain the data on the broadening of the 
distribution at excited rotational levels obtained by ana- 
lyzing the absorption system of QSO 1232+082. Based on 
the calculated populations of rotational levels with depth, 
we assumed an additional mechanism capable of producing 
an effective broadening. This mechanism is related to the 
possibility of describing a molecular cloud by two different 
regions: the shell and the inner part. This separation, which 
naturally arises when the radiative transfer is considered in 
the cloud, stems from the fact that the excited and ground 
rotational levels are mostly populated in the shell and the 
inner part, respectively. In that case, the shell will be a PDR 
in which the influence of ultraviolet photons on the chemi- 
cal composition and thermodynamic balance of the medium 
is strong. In view of the assumptions made in our model, we 
introduced a spatial nonuniformity of the Doppler param- 
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eter b to take into account the PDR. The increase in b in 
the shell is related to the heating of the medium or to the 
stronger turbulence. With a nonuniform parameter b, we 
can obtain the broadenings of the distributions at excited 
rotational levels needed to explain the observational data. A 
similar broadening mechanism for the velocity distributions 
of molecules was used bylNotcrdacm eet al. 20071 However, 
they did not consider the effect that we showed and that is 
related to radiation directionality. We point out that when 
a PDR manifests itself in the observations, the radiation di- 
rectionality effect should also manifest itself. Therefore, in 
the case of a PDR region, it should be taken into account to 
eliminate the systematic error when the observational data 
are analyzed. 



Substituting ([24]) here, we obtain 
ni(v z ,z) = J dv z P ok (v z ,z)x 

— OO 

+00 

x J dot w th (v' z - vt)wturb(vt)w th (v z - v t ) . (26) 

Integration over v t yields the product of twoGaussian 
functions with parameters b and b + = + b 2 h — 

b t h^l + {hurb/b) 2 . 

Equations and ((23} then give ©. 
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APPENDIX 

To derive the distribution thermalization while the turbu- 
lent component is retained, we will assume that that the 
distribution at the ground level in ([5]) is an equilibrium one 
and use the assumption about microturbulence. We can 
then write 



ni(v Zl z) 



n (z ) 
Ai 



Pok{v z ,z) J dv t Wturb{vt)wth(v z - V t ), 

— OO 

(22) 

where no(z) is the number density at the ground level and 
Pok is defined by Eq. (|2|); wturb and wth are the Gaussian 
functions with parameters b tur b and bth , respectively. Here, 
we used the fact that the particle velocity distribution was 
a convolution of the thermal and turbulent distributions. 
Integrating Eq. ([22]) over v z yields 



ni(z) 



n Q {z) 
A, 



-oo -]-oo 



dv t / dv z w turb (v t )w th (v z -v t )P ok (v z ,z). 



— OO — OO 



(23) 

We thus see that the distribution of the upward thrown 
particles in turbulent velocity is 

Wturb(Vt,z) = —--«) fa rt(v t )X 

ni{Z) j±i 

+00 

X J dv zWth (v z -v t )P 0k (v z , 2 ). (24) 

— OO 

Under our assumptions, this nonequilibrium turbulent 
distribution is retained. Thus, after the thermalization of 
the medium, the particle distribution at level I can bewrit- 
ten as 
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